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ABSTRACT 



For the first time, we systematically explored the population of discrete X-ray sources in the outskirt of early-type galaxies. Based 
on a broad sample of 20 galaxies observed with Chandra we revealed overdensity of X-ray sources in their outskirts. They appear as 
halos of resolved sources around galaxies, distributing much broader than the stellar light, and extended out to at least ~ lOr^ (r^ is 
the eff'ective radius). These halos are composed of sources fainter than ~ 5 x 10^^ erg/s, whereas the more luminous sources appear 
to follow the distribution of stellar light, suggesting that the excess source population consists of neutron star binaries. Dividing the 
galaxy sample into four groups according to their stellar mass and specific frequency of globular cluster, we find that the extended 
halos are present in all groups except for the low mass galaxies with low globular cluster content. We propose that the extended halos 
may be comprised of two independent components: (i) LMXBs located in blue (metal-poor) globular clusters (GCs), which GCs are 
known to have a broader distribution than the stellar light; (ii) neutron star LMXBs kicked out of the main body of the parent galaxy 
by the supernova explosion. The available deep optical and X-ray data of NGC 4365 supports this conclusion. For this galaxy we 
identified 60.1 + 10.8 excess sources in the (4 - 10)re region of which ~ 40% are located in globular clusters, whereas ~ 60% are 
field LMXBs. We interpret the latter as kicked NS LMXBs. We discuss implications of these results for the natal kick distributions of 
black holes and neutron stars. 

Key words. X-rays :binaries - globular clusters: general - Galaxy: halo - supernovae: general 



1. Introduction 

LMXBs consist of a compact object - either a black hole (BH) or 
a neutron star (NS) - and a low-mass donor star (< 1 M©), which 
transfers mass via Roche-lobe overflow. The infalling matter is 
heated to X-ray temperatures and releases X-ray luminosities 
in the range of 10^^ - 10^^ erg/s. Due to their luminous na- 
ture, LMXBs add a major contribution to the total X-ray emis- 
sion of galaxies, moreover they determine the X-ray appearance 
of relatively gas-poor earl y-type galaxies (e.g. IGilfanovl 120041: 
iBogdan & Gilfanovl l201lh . The populations of LMXBs have 
been extensively studied in e arly-type galaxies \yith Chandra 
observations (e.g. Irwin et al.|[2003: Gilfanovl l2004l: IZhang et al." 
[201 1: Colbert etal. 2004) and a simple picture of the correla- 
tion of their properties with the parameters of the host galaxy 
emerged. It has been established that LMXB population scales 
with the stellar mass of the host galaxy. For luminous LMXBs 
(Lx > 3 X 10^^ erg/s), it also has been shown that their sur- 
face density closely follows the near-infrared light distribution 
^ilfanov 2004), excep t for the very inner regions of galaxies 
Voss& Gilfanovl I2007h . However, in the course of more than 
ten years of operation of Chandra, a large number of early- 
type galaxies has been observed to a fairly deep sensitivity limit 
(Lx ~ 10^^ erg/s), allowing to probe LMXB populations in un- 
precedented details. 

Such a detailed study of LMXB populations has been per- 
formed for Sombrero galaxy (Ml 04) by Li et al. (2010). Their 
results appear to challenge the empirical picture that the LMXB 
populations closely follow the stellar light distribution at all 



galactocentric radii. As opposed to the expectations, the outskirt 
of Sombrero exhibits a significant X-ray source excess. In the 
halo of this galaxy 101 sources were detected, whereas the ex- 
pected number of cosmic X-ray background (CXB) sources is 
52 ± 1 1, implying a ~ 4Acr excess. Although Li et al I (l2010l) did 
not identify the origins of the excess X-ray sources, they consid- 
ered that either supernova kicked binary systems and/or LMXBs 
associated with GCs could be responsible. 

GCs have a bimodal color distribution (Forbe s et"aDll997l) . 
The red (metal-rich) population follows the stellar light distribu- 
tion, whereas the blue (metal-poor) population has broader dis- 
tribution and traces the more extended dark matter halos (e.g., 
Bassino et al. 200^ . Since few percent of blue GCs host LMXBs 
(Kundu et al. 200i lJordan et al.ll2004 ). it is likely that LMXBs 
associated with blue GCs will have a broader distribution than 
the stellar light. This eflTect could result in an extended popula- 
tion of LMXBs in the outskirts of galaxies. Obviously, the im- 
portance of GC-LMXBs is expected to be more prominent in 
galaxies with higher globular cluster specific frequency (S^). 

Neutron stars receive kicks (so called natal kicks) when they 
are formed in the core collapse supernova explosions. The am- 
plitude of the natal kick velocities can be determined, for ex- 
ample, from proper motion of pulsars, the mean birth speed for 
young pulsars is ^ 400 km/s with the fastest n eutron stars moy - 
ing with the velocity in excess of - 10^ km/s (Hobbs et alj2005|). 
Although only a relativel y small fraction of b inary systems can 
survive large natal kicks (iBrandt & Podsiadlowskilll995h . those 
do survive, can travel to large distances producing a distribution 
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Table 1. The galaxy sample. 
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Notes. (1) - Galaxy distance derived by iTonrv et al.l ([200 Ih. (2) - Scale conversion. (3), (4) and (5) - ^s-hand half-light radius, axis-ratio and 
position angle from the 2MASS Large Galaxy Atlas ('Jarrett et al.''2 003 | ). (6) - Total -band luminosity calculated from t he total apparent ^s- 
band magnitude. (7) - ^s-band mass-to-light ratios derived from Bell & de Jonj (1200 Ih . with B - V colors from RC3 catalog dde Vaucouleurs et al] 
1991). (8) - Total stellar mass calculated from ^s-band total luminosity and the ^s-band mass-to-light ratio. (9) - Globular cluster specific frequen- 
cies derived following Zhans et al. ( 201^. R eferences for the nu mber of observed globular clusters in each galaxy are - ^ Kissler-Patis et al. (1996); 
te^itler et al. (2008); ^Forbes et al. (2001); ^Kissler-Patis et al' (^1997); ^Forbes et al. (1998); ^Harris (1991); ^Rhode & Zepf (2004^; ^iHumphrevI 
(12009): Sikkema et al. (2006); ^Forbes (1996); ^Pens et al. (2008); Dirsch et al. (2005); '^Dirsch (1996); ' Cantiello et al. (2007). (10) - Stellar 
a ge of the galaxy. References - ^^ Terlevich & Forbes (2002); ^^Sanchez-Blazquez et al. (2006); "^ AnnibaH et all ([2007); ^"Thomas et al. (2005D; 
^ iSchweizer & Seitz er (1992). (11) - Total exposure time of the combined Chandra observations. (12) - Limiting source detection sensitivity in 
the 0.5 - 8 keV band in the (0.2 - 3)re {U^J and outer (4 - 10)re (L°^) regions. 



much broader than that of stars. The effect of the natal kicks on 
binaries should be stronger in lower mass galaxies, where kick 
velocities may exceed the escape velocity in the gravitational 
potential of the parent galaxy, and binary systems may leave the 
galaxy or be displaced to very large distance from its main body, 
where binaries were formed. 

In this paper, we aim to comprehensively explore the LMXB 
populations in a large sample of galaxy outskirts with the lim- 
iting sensitivity of ~ 10^^ erg s"^. We will address two major 
points. First, we will investigate whether the existence of ex- 
cess LMXBs is ubiquitous in the outskirts of early-type galax- 
ies. Second, we will study the origin of the excess sources with 
a particular focus on GC-LMXBs and supernova kicked X-ray 
binaries. To achieve our goals, we will use a sample of 20 early- 
type galaxies with deep Chandra observations. Since our galaxy 
sample covers a broad range in both stellar mass and globular 
cluster specific frequency, we will be able to perform a system- 
atic study and address the importance of GC-LMXBs and super- 
nova kicked LMXBs. 

The paper is structured as follows. In Sect. [2] we introduce 
the analyzed sample. In Sect. [3] the X-ray and near-infrared data 
analysis is described. In Sect. |4] and Sect. 5, we discuss the ra- 
dial distributions and the X-ray luminosity distributions of the 
detected X-ray sources. In Sect. Owe describe a case study for 
NGC 4365. Our results are discussed in Sect. [7] and we conclude 
in Sect.[8l 



2. The analyzed sample 



In the present work we explore the LMXB populations in the 
outskirts of early-type (E/SO) galaxies. To achieve this goal and 
obtain statistically significant c onclusions , a bro ad sample of 
galaxies must be analyzed. In Zhang et alj (1201 2h we carefully 
built a sample of 20 galaxies, which also well suits the pro- 
posed goals of the present analysis for the following four rea- 
sons. First, the sample consists of galaxies within the distance 
range of 9.7-27.7 Mpc. Therefore, adequate (~ 5 x 10^^ erg s"^) 
source detection sensitivities can be achieved with moderately 
deep (< 150 ks) Chandra exposures, and the sample galax- 
ies and their outskirts fit well within the ACIS field-of-view. 
Second, the sample includes both GC-rich and GC-poor galaxies 
(5'n= 0.50- 12.38). Third, the sample covers a relatively broad 
range of stellar masses (M^ = (6.1 - 35.6) x 10^^ M©), imply- 
ing vastly diflferent numbers of field LMXBs, as well as diflferent 
dark matter halo masses. Finally, the selected galaxies do not ex- 
hibit ongoing or recent star-formation, hence their stellar content 
is fairly homogeneous and the population of X-ray binaries is not 
polluted by high-mass X-ray binaries. The physical properties of 
the sample galaxies are listed in Table [TJ 
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Table 2. The statistics of X-ray point sources. 
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Notes. (1), (2) and (3) - Number of detected X-ray sources, estimated 
CXB sources, and estimated LMXBs from ^s-band light in the (0.2 
- 3)re region. (4), (5) and (6) - The same quantities in the (4 - 10)re 
region. 



3. Data analysis 

3.1. Chandra data 

The analyzed Chandra observations and t he main step s of th e 
data analysis agree with those outlined in IZhang et alJ (l2Q12h . 
The data was reduced using standard CIAO threads (CIAO ver- 
sion 4.2; CALDB version 4.2.1). To detect point sources, we 
applied the CIAO wavdetect task with parameters taken from 
[Voss & Gilfanovl(l2006ll2007D . Namely, we set the threshold pa- 
rameter (si g thresh) to 10"^, which implies 1 false detection per 
10^ pixels or one ACIS-S CCD. The source detection was per- 
formed in the 0.5 - 8 keV band. To increase the source detec- 
tion sensitivity, we did not exclude high background periods, 
since they are outweighed by the increased exposure time. We 
produced exposure maps in the 0.5-8 keV band, assuming a 
single power-law model with F = 1.7 and Galactic absorption. 
Several galaxies in our sample have been observed in multiple 
Chandra pointings, ther efore we corrected the offsets following 
IVoss & Gilfa^(l2007h using the CIAO reproject .events tool. 
After correcting the offsets, the images were combined and re- 
analyzed. 

To compute the source net counts, we employed circular 
apertures centered on the central coordinates of each source. The 
radius of the aperture was defined as a circle, which includes 
85% of the point spread function (PSF) value. For each source, 
the PSF was determined using the CIAO mkpsf tool. To account 
for the background components, we applied circular regions 
with three times the radius of the source regions, while the over- 
lapping regions of adjacent source s were excluded. The s ource 
net counts were derived following Voss & Gilfanovl (l2007l) . The 
observed net counts were converted to 0.5 - 8 keV band unab- 
sorbed luminosities assuming a power-law spectrum (F = 1.7) 
with Galactic absorption. 



The source detection sensitivity varies throughout the 
Chandra images due to the varying level of the diffuse X-ray 
emission associated with the galaxy, the PSF deterioration at 
large off-axis angles, and the non-uniform exposure of the (com- 
bined) images. For the combined images, we derived the source 
detection sensitivities by inverting the detec tion method, us- 
ing the local PSF, background, and exposure ( Voss & GilfanoyI 
2006). The incompleteness function, K(L), in a certain area was 
calculated by accumulating the sensitivities of the pixels in- 
cluded in the area, weighted by the assumed spatial distribution 
of sources. Thus, K(L) was computed separately for the CXB 
sources with a flat spatial distribution, and for LMXBs assum- 
ing that they follow the ^s-band stellar light. We used these 
functions to calculate th e number of predicted CXB sources 
(Georga kakis et al.ll2008 "l) and LMXBs (Table [2]), and to obtain 
their predicted radial source density profiles (Sect. 14. ll) . The in- 
completeness function was also used to correct apparent X-ray 
luminosity functions (XLFs), and to produce the final XLFs pre- 
sented in Sect. 15.21 

The contribution of CXB sources is estimated based on their 
luminosity function from Georgakakis et al. ( 2008). To convert 
their 0.5-10 keV band log - log 5" distribution to the 0.5 - 8 
keV band, we assumed a power-law model with F = 1.4. Within 
the central regions of the galaxies, the contribution of CXB 
sources is < 10% (Table [2]). Although in the outskirts of the 
sample galaxies the ratio of predicted CXB sources is signifi- 
cantly higher (Table [2]), in most galaxies a statistically signifi- 
cant source excess is detected above the predicted CXB level. 
We stress that the excess sources cannot be attributed to the 
cosmic variance, which can be responsible for variations on the 
~ 10 - 30% level. The accuracy of CXB subtraction is further 
discussed in Sect. 15.11 



3.2. Near-infrared data analysis 

To trace the stellar light of the sample galaxies, we relied on the 
^s-band (2.16 ^m) images of the 2MASS Large Galaxy Atlas 
(iJarrett et al.l 12003). The provided images are background sub- 
tracted, except for NGC 821 and NGC 1052, where we esti- 
mated the background level using nearby regions. To avoid pol- 
lution by bright foreground and background objects, we visually 
identified and removed them. We converted the observed num- 
ber of ^s-band source counts (S) to apparent magnitudes using 
mK = KMAGZP - 2.5 log S, where KMAGZP is the zero point 
magnitude for the ^s-band image. The thus obtained magnitudes 
were used to compute /^s-band luminosities, assuming that the 
absolute ^s-band magnitude of the Sun is Mk,© = 3.39 mag. 

Based on the ^s-band luminosities and the ^s-band mass- 
to-light ratios, we derived the stellar mass of the sample galaxies 
(Table [TJ. The mass-to-ligh t ratios were derive d from the 5 - V 
color indices (RC3 catalog; de Vaucouleurs et al. 1991) and the 
results of galaxy evolution modeling (Bell & de Jong..200L) . 



4. Excess X-ray sources in galaxy outskirts 

4.1. Radial distributions of all galaxies 

To map the spatial distribution of X-ray sources in and around 
the sample galaxies, we built stacked radial source density pro- 
files (Fig. [T]). The profiles were extracted from concentric el- 
lipses, whose shape and orientations were determined by the ^s- 
band photometry of the galaxies (Table [TJ. The contribution of 
CXB sources was subtracted as described in Sect. 13.11 The ob- 
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Fig. 1. Stacked radial source density profiles of the CXB sub- 
tracted LMXBs in the sample galaxies (crosses). The upper and 
lower panels correspond to all the sources and sources brighter 
than 5 X 10^^ erg/s, respectively. The profiles are not corrected for 
incompleteness. The predicted distributions of LMXBs based on 
the ^s-band light are plotted with solid histograms. They take 
into account the source detection incompleteness, as described 
in the Sect. 13.11 and therefore can be directly compared with ob- 
served profile. An excess of LMXBs beyond ~4re is present for 
all sources, but is absent for bright sources. 



served source density profiles are comp ared with those ex pected 
based on the average XLF of LMXBs ( Zhang etal]l2012h . 

In the upper panel of Fig. [T] we show the source density 
profile for all X-ray sources without employing a luminosity 
cut. Within the (0.2 - 3)re region the observed and predicted 
radial profiles are in fairly good agreement with each other. 
Contrarily, in the (4 - 10)re region a highly significant source 
excess is detected. At the first sight , this resu lt appears to con- 
tradict the source density profile of iGilfanovl (120041 who con- 
cluded that the distribution of LMXBs closely follows the stel- 
lar light. However, at the time of writing, only relatively shal- 
low C handra observations of early-type galaxies were available, 
hence IGilfanovl (i2004) considered only bright X-ray sources. In 
the lower panel of Fig. [T] we present the radial source density 
profile including only sources brighter than 5x10^^ erg/s. The 
observed distribution of bright sources follows the stellar light 
distrib ution at all radii, being in good agreement with Gilfanov 
(l2QQ4l) . 

Figure [T] clearly shows that there is an excess in the number 
of faint sources at large radii. The origin and characteristics of 
the faint excess sources are discussed throughout the paper. 

4.2. Source excess in individual galaxies 

In Fig. [T] we demonstrated the existence of an excess population 
of faint sources beyond ~ 4re. However, the presented stacked 
source density profiles were obtained for 20 galaxies, hence it 



reflects the general behavior of our sample. To study individual 
galaxies, we derived the number of detected and predicted X-ray 
sources in each galaxy separately. Given the relatively low num- 
ber of detected sources in individual systems, we did not attempt 
to build radial source density profiles for each galaxy. Instead, 
we divided the galaxies into two regions: the inner region was 
defined as an elliptic annulus with (0.2 - 3)re radius, whereas the 
outer region is an elliptic annulus with (4 - 10)re. Note that the 
central 0.2re region was excluded to avoid the possible source 
confusion. The results are summarized in Table [2l where we list 
the numbers of detected sources, predicted LMXBs, and esti- 
mated CXB sources in the inner and outer regions. As before, 
the predicted numbers of LMXB s were derived using the average 
LMXB XLF (IZhang et al.ll201lh and the stellar mass enclosed in 
the regions. The estimated CXB level was derived as described 
in Sect.O 

In the inner region the numbers of detected and predicted 
X-ray sources agree within ~ 50% for all galaxies, except 
for NGC 4382, NGC 4636, and NGC 4649, where they agree 
within factor of ~ 2. In the entire sample, the total number of 
LMXBs in the inner region after CXB subtraction is 1327.6, 
whereas the expected number is 1217.3, which numbers agree 
within ~ 10%. The total stellar mass in the (0.2 - 3)re region 
is 2.1 X 10^^ Mq, implying a sample-averaged LMXB specific 
freq uency of 6 . 3 + 0. 2 sources per 10^^ M©, in good agreement 
with IGilfanovl (|2004 : IZhang et al.1 (12012 ). In the outer region, 
14 galaxies out of 20 show a notable X-ray source excess. In 
most galaxies the number of additional X-ray sources exceed the 
predicted number of LMXBs by factor of 2.5 - 6, and in NGC 
4636 and NGC 4649 this ratio is even larger. In the outer re- 
gion we detected 609.4 X-ray sources after subtracting the CXB 
level, whereas 158.1 sources are expected, implying a factor of 
~ 3 excess. The total enclosed stellar mass in the outer region is 
2.1 X 10^^ Mq, hence the sample-averaged LMXB specific fre- 
quency is 29.0 ±1.2 sources per 10^^ M©, which is ~ 3.5 times 
more than in the inner region. 

Based on the study of individual galaxies, we obtain two ma- 
jor conclusions. First, the existence of excess X-ray sources ap- 
pears to be a general phenomenon in the outskirts of galaxies. 
Second, the LMXB specific frequency is, on average, by a fac- 
tor of ~ 3.5 higher in the outer region than in the inner region, 
indicating that the excess X-ray sources do not originate from 
CXB fluctuations, and are not directly associated with the stellar 
populations in the outskirts of studied galaxies. 

4.3. The role of supernova kicks and GC-LMXBs 

There are at least two plausible factors which could be respon- 
sible for the observed excess in the surface density of X-ray 
source excess in the galaxy outskirts. They could be LMXBs 
residing in globular clusters and/or accreting NS binaries that 
were kicked to large galactocentric radii by the supernova explo- 
sion, which produced the NS. The importance of GC-LMXBs is 
higher in galaxies with large globular cluster specific frequency 
(5'n), whereas supernova kicked LMXBs are more likely in mas- 
sive galaxies, which have a larger number of potential X-ray 
binaries and can bound the kicked systems due to their deeper 
potential well. To probe whether GC-LMXBs and/or kicked 
LMXBs play a role in the galaxy outskirts, we divided the sam- 
ple into four diff'erent groups, each of them consisting five galax- 
ies. We used 2.0 to separate galaxies with low/high 5'n, and 
= 1.2 X 10^^ Mq to distinguish low/high mass galaxies. The 
thus obtained four groups are: 1) large and large 5'n, 2) large 
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Fig. 2. Stacked radial source density profiles of the CXB subtracted LMXBs in four groups of galaxies with diff'erent and S^, 
marked in the plots. The profiles are not corrected for incompleteness. The predicted distributions of LMXBs based on the ^s-band 
light are plotted with solid histograms. They take into account the source detection incompleteness, as described in the Sect. 13.11 
All groups show excess LMXBs in their outskirts, except for galaxies with small and small S^. 



Mi, and small 5'n, 3) small and large 5'n, and 4) small 
and small S^. 

For each group, we derived stacked radial profiles of the 
detected X-ray sources and computed the predicted number of 
LMXBs as described in Sect. 14.11 The obtained profiles are de- 
picted in Fig. [21 A common property of the four groups is that in 
the inner region the observed and estimated profiles are in fairly 
good agreement. However, in the outer region the first three 
groups show a significant source excess, with the only excep- 
tion being low-mass galaxies with low S^. Thus, we conclude 
that the X-ray source excess cannot be only associated with ei- 
ther GC-LMXBs or supernova kicked LMXBs, but both types of 
sources play a role. 



5. Combined XLFs of LMXBs 

5.1. Accuracy of the CXB subtraction 

To further explore the properties of LMXBs in the sample galax- 
ies, we built combined XLFs of all the sources residing in the 
inner (0.2 - 3)re and outer (4 - 10)re regions. In total, we re- 
solved 1412 and 1077 point sources in the inner and outer re- 
gions (Table [2]), allowing us to build accurate XLFs. 



Since our study is focused on galaxy outskirts with rela- 
tively low source density, it is crucial to accurately account for 
the CXB level. In the inner regions the estimated average con- 
tribution of CXB sources is only ~ 6%, whereas it is ~ 43% 
in the outer regions (Table O. Given that the CXB source den- 
sity exhibits 10 - 30% field-to-field variations due to the cosmic 
variance, the inaccurate subtraction of the CXB level could sig- 
nificantly influence our study. Therefore, in Fig. [3] we show the 
combined cumulative luminosity distributions of the detected X- 
ray sources in the inner and outer regions, which are compared 
with the predicted, incompleteness-corrected CXB source dis- 
tributions. The upper panel of Fig. [3] demonstrates that in the 
inner region the contribution of CXB sources is fairly low at lu- 
minosities below 10^^ erg s"^. Above this threshold we detected 
ultra-luminous X-ray sources, which are most likely stellar mass 
black holes accreting from a low or intermediate mass compan- 
ion. For a comprehensi ve discussion of th ese ultra- luminous X- 
ray sources we refer to lZhang et al.l (l2012h . 

The lower panel of Fig. [3] shows that in the outer region the 
luminosity distribution of sources above 5 x 10^^ erg/s is in ex- 
cellent agreement with the predicted CXB level. Indeed, above 
this luminosity 7 1 ± 8.4 X-ray sources are detected, whereas 65. 1 
CXB sources are predicted. Note that in the bright end the XLF 
is not afl'ected by incompleteness efl'ects. Below 5 x 10^^ erg/s 
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Fig. 3. Cumulative luminosity distributions of all resolved point 
sources in the inner (upper) and outer (lower) regions of the sam- 
ple galaxies. The distributions are not corrected for incomplete- 
ness or the contribution of CXB sources. The shaded areas corre- 
sponds to Icr Poissonian uncertainties. The s olid lines show the 
predic ted CXB luminosity distributions from Georgakaki s et al.l 
(l2Q08h . modified by the incompleteness function, as described in 
the Sect.O 



the number of detected sources is more than a factor of two of 
the expected number of CXB sources. Thus, in agreement with 
the radial source density profiles (Fig. [B, the XLF also shows 
the presence of an excess source population with luminosities 
< 5 X 10^^ erg/s in the outer region. The good agreement between 
the observed source and predicted CXB luminosity distributions 
in the bright end demonstrates the accuracy of our CXB subtrac- 
tion procedure. Moreover, the large diff'erence between the XLFs 
in the faint end excludes the possibility that the observed excess 
is due to the cosmic variance. 

5.2. Source XLFs of the inner and outer regions 

We directly compare the combined XLFs of all galaxies in the 
inner and outer regions in Fig. |4l To build the XLFs, we only 
considered detected X-ray sources above the 0.6 incomplete- 
ness level. The sources observed in individual galaxies were 
combined and weight ed with the enclosed stellar mass within 
the regions following IZhang et aP (l2Qlll) . We subtracted the 
CXB level and applied incompleteness correction as described 
in Sect. 13. 11 In the upper and lower panels of Fig. |4] we show the 
final cumulative and diff'erential forms of the XLFs, respectively. 

Although statistical uncertainties at high luminosities are 
rather large, the XLFs of the inner and outer regions appear to 
be consistent with each other in the bright end. This conclusion 
is further supported by the luminosity distributions shown in the 
Fig. [3] and the shape of the radial profile of bright sources in the 
Fig. [TJ On the other hand, the outer region shows a significant 
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Fig. 4. XLFs of LMXBs in the inner and the outer regions in cu- 
mulative (upper panel) and diff'erential (lower panel) forms. The 
data for the outer regions is marked by circles in the lower panel 
and is surrounded by the shaded area showing the Icr Poissonian 
uncertainty in the upper panel. The statistical uncertainties for 
the inner regions have a smaller amplitude. 



source excess in the luminosity range below ~ 5 x 10^^ erg/s. 
This threshold luminosity is fairly close to the Eddington lumi- 
nosity of an accreting NS. Therefore, it is reasonable to assume 
that the bulk of faint excess sources are NS binaries. This con- 
clusion is consistent with th e proposed origin of excess sources . 
Indeed, both GC-LMXBs (Portesies Zwart & McMillan "20001) 
and supernova kicked binaries (Brandt & Podsiadlowski 1995) 
are predominantly NS binaries. 

In principle, it would be beneficial to compare the XLFs of 
the four groups of galaxies with low/high and S^. The shape 
of the XLFs could help to constrain the origin of the excess 
sources, since field and GC-LMXBs exhibit markedly diff'erent 
XLFs at luminosities below few times 10^^ erg/s (Zhang et all 
|201lV However, the studied sample is not suitable for such a 
comparison since the source detection sensitivity of most sam- 
ple galaxies is comparable to this luminosity limit. Above few 
times 10^'^ erg/s the XLFs of field and GC-LMXBs are virtually 
indistinguishable, hence we do not attempt a more detailed study 
of luminosity distributions. 

6. A case study for NGC 4365 

The stacked source density profiles (Fig. [2]) indicate that the ob- 
served X-ray source excess in the galaxy outskirts cannot be ex- 
plained only by the population of GC-LMXBs, but most likely 
supernova kicked accreting NSs also play a notable role. To 
separate the population of GC-LMXBs and supernova kicked 
sources, GC-LMXBs must be identified at large galactocentric 
radii. However, for most galaxies in our sample this is not fea- 
sible due to the limitations of the available globular cluster data. 
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Although several galaxies in our sample have been observed by 
the Hubble Space Telescope (HST), these images do not cover 
their outskirts. Ground based observatories have larger field-of- 
view, but these images usually suffer from major incompleteness 
and foreground star contamination. An exception from this trend 
is NGC 4365, for which deep optical data is available, allowing 
to detect GCs even in its outskirts. Additionally, deep Chandra 
data is also available with a limiting luminosity of ~ 10^^ erg/s. 
Therefore, we probe the X-ray populations of NGC 4365 as a 
case study. 

NGC 4365 is a large elliptical galaxy with M^= 1.8 x 
10^^ M© and S^= 3.95, which exhibits a notable source excess 
in the (4 - I0)r^ region. The GC population of NGC 4365 has 
been studied as part of the SAGES Legacy Unifying Globulars 
and Galaxies Survey (SLUGGS). The galaxy and its outskirts 
have been observed in 8 pointings with the Advanced Camera 
for Surveys (ACS) onboard HST, which data has been comple- 
mented with ^t /Z^am/S-Cam obser vations. Based on this com- 
bined data set, iBlom et aP (l2Q12h reported more than -6000 
GC candidates in and around NGC 4365. Beyond the central 
0.5' (~lre) the obtained list of GC candidates is complete to 
the turnover magnit ude, additionall y not significantly aff'ected 
by contamination (Blom et al.ll2012h . 

To build the list of GC candidates around NGC 4365, we 
relied on the publicly available SLUGGSQ data set, combining 
the HST and Subaru identifications together. If a GC candidate 
has been observed by both telescopes, we prioritized HST. We 
identified GC-LMXBs by cross-correlating the coordinates of 
GC candidates with the coordinate s of detected X-ra y sources 
following the method described in Zhang et aP (12011 ). The ap- 
plied match radius was 0.5'', which results in a total number of 
104 X-ray sources residing in GCs within the study field. The 
estimated number of random matches is ~ 3. The 104 identified 
sources are designated as GC-LMXBs, whereas all the others are 
field sources. 

In Fig. [5] we separately show the radial source density pro- 
files of GC-LMXBs and CXB subtracted field sources, which 
are compared with the expected number of LMXBs based on the 
enclosed stellar mass and the average LMXB XLF (Sect. 14.21) . In 
the (4 - 10)re region 116 X-ray sources are detected (Table O, 
which significantly exceed both the CXB prediction (44.6) and 
the expected number of LMXBs associated with the stellar light 
(1 1.3). Among the 116 sources, 25 sources are identified as GC- 
LMXBs, implying that 46.4 ± 9.5 field X-ray sources are associ- 
ated with the galaxy, which is ~ 3 times higher than the predicted 
LMXB number. Thus, we estimate that 35.1 ± 9.5 sources could 
be associated with supernova kicked LMXBs. 

Based on the above interpretation we conclude that 1) the 
source excess in the outskirts of NGC 4365 cannot be only at- 
tributed to GC-LMXBs, but presumably supernova kicked NS 
binaries also contribute; and 2) the importance of supernova 
kicked NS binaries is comparable to that of GC-LMXBs in a 
massive galaxy even with relatively high S^. 



7. Discussion 

7.1. LMXBs associated with intracluster light 

In galaxy clusters the stellar light is not only associated with 
luminous galaxies, but certain fraction (~ 5 - 50%) of the to- 
tal optical cluster luminosity is associated with the intraclus- 
ter light (ICL), which is bound to the cluster potential (e.g.. 
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Fig. 5. Stacked radial profiles of sources in NGC 4365. The 
CXB subtracted field LMXBs are marked with squares. The GC- 
LMXBs are marked with circles. The predicted distribution of 
LMXBs from ^s-band light are plotted with histograms. 
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iFeldmeier et ani2004l: iGonzalez et aDl2005l: IZibetti et al.ll2005l) . 
Due to its faint surface density , the ICL can onl y be detected 
with deep optical observations ( Mihos et"ani2005l) . and remains 
essentially undetected in 2MASS images. Similarly to the stellar 
bodies of galaxies, LMXBs can evolve in the ICL, which in prin- 
ciple could contribute to the observed source excess in the galaxy 
outskirts. Since part of the studied galaxies reside in rich galaxy 
cluster environment, we investigate the importance of LMXBs 
associated with the ICL. 

To estimate the average surf a ce de nsity of the ICL, we relied 
on the analysis of IZibetti et al.l (1200 5), who stacked the Sloan 
Digital Sky Survey imaging data of 683 clusters in the redshift 
range of z = 0.2 - 0.3. We used their stacked /-band ICL sur- 
face brightness distribution and assumed an /-band mass-to-light 
ratio of M^/Lj = 1.50 to deduce the stellar mass. The aver- 
age ICL mass density profile extending to 1000 kpc is shown 
in the upper panel of Fig. [6j To convert the stellar mass den- 
sity profiles to the predicted LMXB source density profile, we 
assumed a galaxy cluster at a distance of 15 Mpc, adequate for 
Virgo cluster, and a Chandra exposure time of 100 ks, which 
corresponds to a detection sensitivity of ~ 10^^ erg/s. Using 
these parameters and the average ICL mass profile, we com- 
puted the num ber densities of the predicted LMXBs using their 
average XLF (IZhang et al.lf2012) . The obtained source density 
profile reveals an LMXB density of ~ 10"^ - 10"^ arcmin"^ 
in the radial range of 100 - 1000 kpc (lower panel of Fig. O. 
Note that several galaxies in our sample lie at even larger radii 
from the cluster center, implying an even lower LMXB density. 
Assuming the same distance and Chandra exposure time, we 
estimated the CXB source dens ity from their average log(N)- 
log(S) (Georgaka kis et al.ll2008l) . We obtained an average level 
of ~ 1 arcmin"^, which is 2 - 3 orders of magnitude higher than 
the expected LMXB number associated with the ICL. Since the 
expected LMXB density from the ICL is only a minor fraction 
of the CXB level, and the CXB level exhibits 10 - 30% field-to- 
field variations, intracluster LMXBs are nearly impossible to be 
detected and identified. Thus, LMXBs associated with the ICL 
will not add a notable contribution to the source excess detected 
in the galaxy outskirts. 
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Fig. 6. Up per panel: stel l ar ma ss density profile of the ICL calcu- 
lated from Zibet ti et alJ (|2005). Lower panel: estimated number 
density of LMXB s associated with the ICL (triangles) and CXB 
sources (solid line), assuming a cluster at a distance of 15 Mpc 
and a detection sensitivity of 10^^ erg/s. 




Fig. 7. Stacked radial profiles of the CXB subtracted LMXBs in 
the 10 younger (upper panel) and 10 older (lower panel) galax- 
ies. The estimated numbers of LMXBs from ^s-band light are 
plotted with solid histograms. The excess LMXB populations in 
the outer regions are present for both younger and older galaxies. 



7.2. Stellar age dependence 

In lZhang et al.l (|2012 ^ we demonstrated that a correlation exists 
between the LMXB population and the stellar age of the host 
galaxy. Namely, the specific frequency of LMXBs is ~ 50% 
higher in older galaxies than in younger ones. We demonstrated 
that this excess could not be explained by only the fact that older 
galaxies have larger GC-LMXBs populations, but also intrinsic 
evolution of the LMXB population with time is likely to play a 
role. In the view of these findings we investigate whether the ob- 
served source excess in the outskirts of galaxies correlates with 

the stellar age^ 

Folio winglZhang eLal.l (|20J2|), the sample of 20 galaxies was 
divided to younger and older galaxies using median age of 6 
Gyrs, which resulted in 10 younger and 10 older systems. We 
built stacked radial source density profiles for both sub- samples 
as described in Sect. 14. II and depicted them in Fig. [71 Compared 
to the predicted LMXB level, both young and old galaxies have 
a significant excess of X-ray sources in the (4 - 10)re region. 
In young galaxies 199.8 sources are observed above the CXB 
level in the outer region, which is twice larger than the predicted 
number of LMXBs from the Ks-hand fight (66.8). For the old 
sub-sample we detect 409.6 sources above the CXB level, being 
a factor of 3.5 times more than the ^s-band prediction (91.3). 
Thus, the source excess appears to be somewhat more signifi- 
cant in old galaxies than in young ones. This result can be in- 
terpreted as the result of two combined eff'ects. First, a strong 
correlation exist s between the stel lar age and the specific fre- 
quency of GCs dZhang et al.l 120121) . Indeed, in our sample the 
median is significantly higher for old galaxies (S^= 4.64) 
than for younger ones (S^= 1.56), implying that the frequency 
of GC-LMXBs is also notably higher for older galaxies in the 
outer regions. Second, if younger galaxies indeed have smaller 



populations of primordially formed LMXBs, there would also be 
fewer kicked sources in the outskirts. Note that the median stellar 
mass of the young and old sample is comparable (1.53x10^^ M© 
and 1.47 xlO^^ M©, respectively), hence the importance of su- 
pernova kicks should be comparable. 

However, due to the correlation between the stellar age and 
5'n, and the lack of complete GC catalogs (Sect. [6]), the popula- 
tion of GC-LMXBs cannot be directly separated from supernova 
kicked neutron stars. Therefore, based on the present sample, it 
is not feasible to comprehensively study the frequency of super- 
nova kicked sources as a function of the stellar age. 

7.3. Quantifying the nunnber of excess sources 

To statistically probe the origin of excess X-ray sources in the (4 
- 10)re region, we characterized the LMXB specific frequency 
in each galaxy using: 

A^x-A^cxB 

/XLF ~ \^) 

M^xf F{L)KiMX^{L)dL 

where Nx and A^cxb are the numbers of resolved X-ray sources 
and the predicted CXB sources. For each galaxy, these values 
agree with those listed for the outer regions in Table [2l is 
the stellar mass in units of 10^^ M©, F{L) is the average diff'er- 
ential XLF normalized to 10^^ M©, and ^lmxb(^) is the incom- 
pleteness function for LMXBs in each galaxy. Note that this for- 
mula is essentially identical with that introduced by lZhang et al.l 
(2012). 

According to equation ([T]), if the number of X-ray sources 
agrees with that predicted from the average XLF, we expect to 
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observe /xlf = 1- However, for the studied 20 galaxies, we ob- 
tained the median value of /xlf = 3.31 in the outer region, in- 
dicating a significant source excess. Assuming that the X-ray 
source excess originates from supernova kicked binaries and 
GC-LMXBs, we fit the data with a two parameter linear model: 

/xlf = ^xM^ -h/?x^N (2) 

where is the stellar mass in units of 10^^ M© and is the 
globular cluster specific frequency. To find the best-fit values, 
we used;^^ minimization and obtained a = 0.081 ± 0.025 and 
b = 0.60 ± 0.12. The large x^/d.o.f = 104.9/18 implies no- 
table scatter in the relation. Taking into account the median stel- 
lar mass (M* = 1.3 x 10^^ M©) and the median specific GC 
frequency (S^ = 2.0) of our sample, we find that (on average) 
the contribution of kicked binaries (a x M^,) and GC-LMXBs 
(b X S^) is comparable. Thus, our statistical analysis confirms, 
that the excess X-ray sources have two diff'erent origins, namely 
from GC-LMXBs and supernova kicked binaries. 

7.4. The role of supernova kicks 

The fate of the kicked binary depends on the amplitude of the 
kick velocity with respect to the orbital velocity in the binary 
system (Vorb) and the escape velocity in the gravitational poten- 
tial of the parent galaxy (Vesc). If the final velocity attained by 
the binary system (Vsys) exceeds the escape velocity, the binary 
will become unbound to its parent galaxy. On the other hand, if 
the kick velocity is too large with respect to the orbital veloc- 
ity in the binary system, it will be destroyed by the kick (e.g. 
iBrandt & Podsiadlowski 1995). Due to this selection eff'ect, the 
average system velocities of LMXB s are expected to be smaller 
tha n for isolated neutron star s . 

iBrandt & PodsiadlowsO (Il995h studied the eff'ect of the 
supernova kick on the binary system. Assuming that the 
birth velocities of neu tron stars are distributed according to 
iLyne & Lorime j (1 19941) with the mean value of 450 km/s, they 
concluded that LMXBs that remain bound after the supernova 
explosion have an average system velocity of Vsys = 180 ± 
80 km s"^. The typical escape velocities of o ur sample galaxie s 
are in the range of Vesc ~ 250 - 1000 km s"^ (IScott et al.ll2009h . 
larger than the average system velocity, therefore we should not 
expect that a significant fraction of LMXBs escape the parent 
galaxy, especially the more massive ones. Nevertheless, sys- 
tem velocities are comparable to the escape velocities, suggest- 
ing that natal kicks of the compact objects can significantly 
modify the spatial distribution of LMXBs in the galaxy, lead- 
ing to its br oadening as compared to the distribution of stars 
(iBrandt & Podsiadlowsk i 1995). However, this does not explain 
the observed dependence on the stellar mass. Indeed, assuming 
that the stellar mass is a proxy to the total gravitating mass, we 
should expect that the eff'ect of kicks is stronger in the lower 
mass galaxies, opposite to what is observed. 

However, it has been suggested that the birth velocity distri- 
bution of Galactic radio pulsar population has a bimodal struc- 
ture, with the low velocity peak at ^ 90 km/s and the high ve- 
locity peak at 500 km/s and the population roughly equally di- 
vided between the two components ( Arzoumanian e t al .112002a) . 
The bimodal distribution of the natal kicks will result in 
the similarly bimo dal velocity distribution of LMXBs (e.g., 
iRepetto et al.ll2012l) . Thus, one may expect that the lower ve- 
locity component remains bound in galaxies of any mass from 
our sample. The fate of the high velocity component, however, 
depends on the gravitational mass of the galaxy. In the higher 



mass galaxies it will remain bound but would be (much) broader 
distributed than stars, due to much larger average velocity, thus 
giving rise to the excess of X-ray sources in the galactic out- 
skirts. In the case of the lower mass galaxies, with the escape 
velocity ^ 500 km/s, the majority of the high velocity systems 
will escape the galaxy and no extended halo of LMXBs will ap- 
pear. 

The escaping binaries will travel to large distances in a rel- 
atively short time scale. Assuming a system velocity of Vsys = 
200 km s"^ and a radial trajectory, we estimate that in 1 Gyr the 
binary will move from the host galaxy to a distance of ~ 200 
kpc. This value significantly exceeds the field-of-view of the 
presently studied Chandra observations (~ 30 - 60 kpc), im- 
plying that such kicked X-ray sources cannot be detected in our 
sample. On another note, due to the expected low surface density 
of supernova kicked sources at large radii, they will be virtually 
indistinguishable from the CXB level. 

As the exact cause of the natal kick is not understood, it is not 
known what birth velocities of black holes are. Unlike neutron 
stars, they cannot be measured directly from proper motions of 
isolated black holes, while attempts to derive them from the spa- 
tial distribution of black hole binaries so far produced ambigu- 
ous results. Among diff'erent models, a momentum-conserving 
model is often considered, which assumes that during the su- 
pernova explosion, black holes and neutron stars receive same 
momentum (rather than same kick velocity). In this model, the 
kick velocities of black holes are by a factor of ~ 5 - 10 smaller 
that those received by the neutron stars, in accord with their mass 
ratios. Therefore black holes kicks would be insufficient to drive 
them out of even the lowest mass galaxies in our sample. This 
agrees with the fact, that the threshold luminosity above which 
the extended halo of X-ray sources vanishes, ~ (3 - 5) • 10^^ 
erg/s, is close to the Eddington luminosity limit for the neutron 
star. It also suggests that the excess sources in the galactic out- 
skirts with supernova kicked origins are neutron star binaries. 

Above, we outlined a qualitative picture which may become 
a framework for a quantitative analysis. In order for this pic- 
ture to be proven or falsified, realistic calculations need to be 
done, similar to the ones performed by Brandt & Podsiadlowski 
d 19951) : iRepetto et al. (2012). Results of such calculations can 
be compared with the observed radial source density profiles in 
galaxies of diff'erent masses. Such a comparison would also re- 
quire detailed globular cluster information in order to remove 
contribution of globular cluster LMXBs. This, however, is be- 
yond the scope of this paper. 

8. Conclusion 

In this paper, for the first time, we systematically explored the 
population of LMXBs in the outskirts of early-type galaxies. 
We studied a sample of 20 early-type galaxies based on archival 
Chandra data, which allowed us to perform a statistically signif- 
icant study of the X-ray populations with a limiting luminosity 
of ~ 10^^ erg s"^ . Our results can be summarized as follows. 

1. We demonstrated the existence of an excess X-ray source 
population in the outskirts of early-type galaxies. These 
sources form a halo of compact X-ray sources around galax- 
ies extending out to at least ~ lOr^ and probably much fur- 
ther. Their radial distribution is much broader that the distri- 
bution of the stellar light. The extended halo is comprised of 
sources with luminosity ^ 5-10^^ erg/s, while the more lumi- 
nous sources appear to follow the distribution of the stellar 
light without any notable excess at large radii. This suggests 
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that the majority of the excess sources are neutron star bina- 
ries. 

2. Dividing the sample galaxies into four different groups based 
on their stellar mass and specific frequency of globular clus- 
ters, we found that the extended halo of compact sources is 
present in all groups except for the low-mass galaxies with 
low globular cluster content (small S^). 

3. We performed a case study of NGC 4365, for which galaxy 
deep optical (HST and Subaru) and X-ray data (Chandra) 
is available, allowing the identification of GC-LMXBs. In 
the (4 - 10)re region we detected 60.1 ± 10.8 excess X-ray 
sources, out of which 25 are GC-LMXBs and 35.1 ± 9.5 are 
field sources which are presumably supernova kicked NS bi- 
naries. 

4. Interpreting our findings, we proposed that the excess 
sources are comprised of two independent components: (i) 
dynamically formed sources in blue (metal poor) globular 
clusters, which are known to be distributed broader than the 
stellar light; (ii) primordial (field) neutron star LMXBs ex- 
pelled from the main body of the galaxy due to kicks re- 
ceived during the supernova explosion. 
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